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Abstract 

We report observations of the Cabibbo-suppressed decays B — > D^K~ 
using a 10.4 fb" 1 data sample accumulated at the T(45) resonance with 
the Belle detector at the KEKB e + e~ storage ring. The high-momentum 
particle identification system of Belle is used to isolate signals for B — > D°K~ , 
D + K~, D*°K~ and D* + K~ from the B — > D^*\~ decay processes which 
have much larger branching fractions. We report ratios of Cabibbo-suppressed 
to Cabibbo-favored branching fractions of: 

B(B~ -» D K~)/B(B- -» D vr-) = 0.079 ± 0.009 ± 0.006; 

,8(5° -» D+K-)/B(B° -> D+7T-) = 0.068 ± 0.015 ± 0.007; 

-» D*°K-)/B(B- -» D*°7r-) = 0.078 ± 0.019 ± 0.009; and 

£(i? -» D* + K~)/B(B° -» D*+7r-) = 0.074 ± 0.015 ± 0.006. 

The first error is statistical and the second is systematic. These are the 
first reported observations of the B — > D + K~, D*°K~ and D* + K~ decay 
processes. 

PACS numbers: 12.15.Hh, 13.25.Hw 
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Comprehensive tests of the Standard Model mechanism for CP violation will ultimately 
require measurements of the 3 angle of the Kobayashi-Maskawa unitarity triangle For 
this, many authors have proposed the measurement of direct CP- violating asymmetries 
due to the interference between b — > c and b — > u transition amplitudes in the Cabibbo 
suppressed B~ — > D°K~ channel as a theoretically clean way to determine 03 0. As a first 
step in this program, it is necessary to establish that the Cabibbo suppressed decay modes 
exist and occur at the expected rates. 

In the tree-level approximation, the branching fractions for the Cabibbo suppressed decay 
processes B — > D^K~ are related to those for their B — > D^n~ counterparts []3| by 

K(n n(*) x~\ 

R - l\ B - - tan 2 9a{f K /Uf - 0.074. (1) 

Here 9q is the Cabibbo angle, and fx and f n are the meson decay constants. The 
only Cabibbo suppressed B — > DK decay observed to date is B~ — > D°K~, 
which is reported by the CLEO group to have a ratio of branching fractions R = 
B(B~ -> D°K-)/B(B- -> D°ti-) = 0.055 ± 0.014 ± 0.005 §, in agreement with expec- 
tations. 

In this paper, we report the first observations of the Cabibbo suppressed processes B — > 
D*°K~, B° —> D* + K~ and B° — > D + K~, and a new measurement of B~ —> D°K~, using 
data collected at the T(45) resonance with the Belle detector || at the KEKB asymmetric 
energy e + e~ collider ||. The good high momentum particle identification capability of the 
Belle detector enables us to extract signals for B — > D^K~ that are well separated from 
the more abundant, Cabibbo favored B — > D^n^ processes. The results are based on a 
10.4 fb _1 data sample that contains 11.1 million BB pairs. 

Belle is a general-purpose detector which includes a 1.5 T superconducting solenoid 
magnet that surrounds the KEKB beam crossing point. Charged particle tracking covering 
approximately 90% of the total center of mass (cm) solid angle is provided by a Silicon 
Vertex Detector (SVD), consisting of three nearly cylindrical layers of double-sided silicon 
strip detectors, and a 50-layer Central Drift Chamber (CDC). Particle identification is ac- 
complished by combining the responses from a Silica Aerogel Cerenkov Counter (ACC) and 
a Time of Flight Counter system (TOF) with specific ionization (dE/dx) measurements in 
the CDC. The combined response of the three systems provides K/tt separation of at least 
'2.5a equivalent' for laboratory momenta up to 3.5 GeV/c. For distinguishing the prompt 
kaons and pions in B — >• D^*'h~ (h~ = K~ or 7r~) decays, only the ACC and dE/dx are used 
since the TOF system provides no significant kaon and pion separation at momenta relevant 
to this analysis. A Csl Electromagnetic Calorimeter (ECL) located inside the solenoid coil 
is used for 7/71- detection. 

The B — > D ( -* S) K~ processes have kinematic properties nearly identical to those of B — > 
D^tt^ decays. While the latter processes produce the most significant backgrounds, they 
also provide control samples that we use to establish cuts on kinematic variables, determine 
the experimental resolutions, evaluate the systematic errors, and normalize the results. 

In this analysis, we require, except for the Kg — > 7r + 7r~ decay daughters, that the charged 
tracks have a point of closest approach to the interaction point within ±5 mm in the direction 
perpendicular and ±5 cm in the direction parallel to the beam axis. For each charged track, 
the particle identification system is used to determine a K/tt likelihood ratio P(K/tt) that 
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ranges between (likely to be n) and 1 (likely to be K). We form candidate D° mesons using 
the K~tt + , K~7[ + it and K~tt + 7t + tt~ decay modes, and D + mesons using K~tt + tt + , Ks^ + , 
i^57T + 7r + 7r _ and K~K + tt + decays. (The inclusion of charge conjugate states is implied 
throughout this report.) For the assignment of charged kaons from D decays, we require 
P(K/tt) > 0.3 in reconstructing D° — * K~7t + and all other D° decay modes associated with 
D* + — > D°7T + , otherwise we require P(K/it) > 0.7. Candidate 7r° mesons are reconstructed 
from pairs of 7's, each with energy greater than 30 MeV, that have an invariant mass in 
the range of ±2a (a = 5.3 MeV/c 2 ) of the measured ir° mass value and magnitude of the 
total three-momentum greater than 200 MeV/c. For the slow tt° from the D* — ► Dtt° decay 
we only require that the invariant mass is in the range between —14.4 and +10.8 MeV of 
the 7T° mass. The K$ — > 7r + 7r~ candidates are reconstructed from two oppositely charged 
tracks that form an invariant mass within ±3<r (a = 4.6 MeV/c 2 ) of the measured Ks mass 
value with a vertex which is displaced from the interaction point in the direction of the 
K s momentum. The selected 7r° and K s candidates are kinematically constrained to the 
nominal mass values. 

For each D-meson decay topology, we select particle combinations that have a recon- 
structed invariant mass within ±2.5<7b of the measured D mass value, where o~d is the D 
mass resolution, which varies from 5 to 13 MeV/c 2 depending on the decay mode. After 
selection, the D candidates are subjected to a mass constrained kinematic fit. For all modes 
except for D + — ► Kg7T + , we further reduce backgrounds by a selection on the quality of the 
vertex fit. 

For D*° and D*+ candidates, we use the D*° -> D°ir°, and D*+ -> D°ir+ and D+ir° 
decay channels. We select events where the mass difference between the Dir and D particle 
combinations is within ±3er of the expected value for D* + — > D°tt + , and ±2cr for D*° — > 
D°7T and D* + — > D + tt°. A kinematic fit that constrains the D* mass to its nominal value 
is applied to the events that satisfy the selection criteria. 

When we isolate B — > D^h~ candidates, we use a quantity we call the lab constrained 
mass, Mic, which is the D^h~ invariant mass calculated with the assumption of two equal- 
mass particles from laboratory momenta: M; c = J (Eg b ) 2 — (pb) 2 , where pb is the B candi- 
date's laboratory momentum and E% = -g—(s/2 + P ee • Pb), where s is square of the center 
of mass (cm) energy, Pb is the laboratory momentum vector of the B meson candidate, and 
P ee and E ee are the laboratory momentum and energy of the e + e~ system, respectively. To 
identify B -> D^K~/D^tt~ samples we use the cm energy difference, which is defined as 
AE = -E^>) + E™ — E£™ m , where E^ is the cm energy of the candidate and Ef^ is 
the cm energy of the prompt h~ track calculated assuming the pion mass and E£™ m is the 
cm beam energy. With this pion mass assumption, B — > D^ir - events peak at AE = 0, 
while the D^K~ peak is shifted to AE = —49 MeV. Typical AE resolutions obtained 
for B -> D«7r- and B -> D^K~ are 16 MeV and 19 MeV, respectively, hence we can 
distinguish these two processess by the AE distribution. For further analysis we accept B 
candidates with 5.27 < M lc < 5.29 GeV/c 2 and -0.2 < AE < 0.2 GeV. 

In the case of multiple entries from one event, we choose the one with the smallest value 
of a x 2 th & t is determined from the differences between measured and nominal values of 
Md, Mic and, when appropriate, Mp* — Md and M^o. For the latter, we only consider 
the 7T° from D*° — > D°7i° and D* + — > D + tt° decays. Background events from e + e~ — > qq 
continuum processes are reduced using the normalized second Fox- Wolfram moment 0, R2, 
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and the angle between the sphericity axis of the B candidate and the sphericity axis of the 
rest of the particles in the event, 6 sp h- For the events with D° — > K~it + decays and all 
modes with D* + — > D 0, k + , continuum backgrounds are small and we only require R2 < 0.5. 
This cut retains 96% of the signal and rejects 25% of the continuum. For all other modes, 
we impose the additional requirement of | cos8 sp h\ < 0.75, which retains 75% of the signal 
and rejects 80% of the continuum events. 

We extract B — > D^K~ enriched samples by applying a stringent particle identification 
condition on the prompt h~, namely P(K/tt) > 0.8; B — > D^n^ samples are selected with 
a criterion P(K/tt) < 0.8. The AE distributions for the B — > D7r~ [B — > DK~ enriched] 
samples are shown in Figs. 0(a) and (b) [(c) and (d)]. The corresponding distributions for 
the D* channels are shown in Figs. |2|(a) through (d). 

Kaon and pion identification efficiencies s(K) and e(ir) are experimentally determined 
from a kinematically selected sample of high momentum D* + — > D°tt + , D° — > K~ir + decays 
where the K~ and 7r + mesons from D° candidates are in the same momentum and angular 
region as the prompt h~ particle in B — > D^*'h~ decay (2.1 < p cm < 2.5 GeV/c). With our 
P(K/tt) cut value of 0.8, the efficiencies are e(K) = 0.765 ±0.006 and e(rr) = 0.980 ±0.003, 
and the tc —> K fake rate is 2.0 ± 0.3%. 

In Figs. |T|(c) and (d), and Q(c) and (d), peaks around AE = —49 MeV correspond to 
B — > D^K" decays while peaks around AE = are due to feed- across from misidentified 
D^n~ decays. The area of the feed-across peak is 2.0% of the D^*>ti~ signal yield (in 
Figs. ||(a) and (b), and 0(a) and (b)), which is consistent with the 7r — > K fake rate. 

We determine the numbers of D(*'ir~ events and the line shape parameters by fitting 
the AE distributions for the D^n~ samples of Figs. |I](a) and (b), and |^(a) and (b). We 
represent the signal distributions using two Gaussian functions with different central values 
and widths. The background has two components. One is due to contributions from Dp~ , 
D*p~ , and other 5-meson decay modes |§, which produce the complicated structures mostly 
seen at negative value of AE, and the other is due to continuum events, which populate 
the entire AE region. The shapes of the BB backgrounds are determined using Monte 
Carlo (MC) simulated events ||; those for the continuum backgrounds are taken from the 
AE distributions for events in the M[ c side band regions (5.20 < M[ c < 5.26 GeV/c 2 ). In 
the fit to the D^ir^ sample AE distributions, we allow the peak positions, widths and 
normalization of the signal functions to vary, as well as the normalization of both the BB 
and continuum background contributions. The fit results are shown in Figs. |I|(a) and (b), 
and ^](a) and (b) as solid curves. The resulting numbers of D^ir~ events are given in Table |. 

In the fits to the D^*'K~ enriched AE distributions, the parameters of the two Gaussians 
for the D(*'K~ signal are fixed at the values determined from fits to the D^ti~ samples 
with a kaon mass hypothesis applied to the prompt pion. The relative position of the signal 
with respect to the original D^n^ signal is reversed. This procedure accounts for the 
kinematic shifts and smearing of the AE peaks caused by the incorrect mass assignment. 
For the feed-across from the D^n^ peak, we fix the parameters at the values determined 
from the D^rr~ fits. In these fits, the areas of the D^*'K~ signal and functions for the 
D^n~ feed-across and the scaling factors of the continuum background shapes are allowed 
to vary. The BB background in the D^K~ enriched sample has two components: modes 
which also contribute to the D^n~ sample with one track misidentified as a high momentum 
kaon and other Cabibbo suppressed modes. The normalization of the feed-across from the 
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first component is fixed to the fit result for the BB background in the D^n~ sample 
multiplied by the measured pion fake rate. The contributions from other D^*>K^*> modes 
are determined from a Monte-Carlo simulation assuming that the branching fractions of the 
suppressed modes relative to the corresponding D^*'tt~ /D^*> p" modes are reduced by the 
usual Cabibbo factor. The fit results are shown as solid curves in Figs. |l](c) and (d), and 0(c) 
and (d). The numbers of events in the D^*'K~ signal and D^*\~ feed-across peaks are listed 
in Table |. Also listed in Table | are the statistical significance values for the D^K~ signals, 
which are determined from the likelihood values of fits made to the D^K~ enriched AE 
distributions with the signal yield fixed to zero. For each channel, the statistical significance 
of the signal corresponds to at least five standard deviations. 
Experimentally, the ratio of branching fractions is given by 

N(B -> D&K-) rj(D^ir-) e(n) 
~ N(B -»• DWtt-) X 7](DWK-) X e(K) 

where N and rj denote the numbers of events and detection efficiencies for the indicated 
processes, and e's are the prompt pion and kaon identification efficiencies, respectively. The 
signal detection efficiencies are determined using MC, and t](D^K~) are approximately 
5% lower than r)(D^ir~) due to the decay-in-flight effect of kaons. Particle identification 
efficiencies are determined as described before. 

Since the kinematics of the B — > D^*'K~ and B — > D(*>ir~ processes are quite similar, 
most of the systematic effects cancel in the ratios of branching fractions. The main sources 
of systematic error that do not cancel are the uncertainties in Kj tt identification efficiencies 
and the shape of the background in the AE distributions. Using MC simulations, we esti- 
mate the systematic error of the K identification due to differences in the angle-momentum 
correlations of the D* + and signal samples to be 2%. To estimate the systematic error due 
to the parameterization of the AE distributions, we use several different methods of fitting. 
These include using linear background functions, fixing the BB backgrounds to MC calcu- 
lated values, and forcing the parameters of D^K~ and/or D^*'tt~ signals to deviate from 
their optimal values by ±lcr. The quadratic sums of those values are used as measures of the 
systematic errors from this source. Compared to these, which range from 7.5 to 10.8%, other 
sources of systematic errors are determined to be negligibly small. The total systematic error 
for each channel is taken to be the sum in quadrature of the individual components. 

The resulting R ratio measurements are listed with their statistical and systematic errors 
in Table |I[ For the B — > D + K~, D*°K~ and D* + K~ decay processes, these are first 
measurements. In all cases, the R ratios are consistent with the expected value given in 
Eq. (1). 

In summary, by taking advantage of the good high momentum particle identification 
capability of Belle, we have observed signals for the Cabibbo suppressed decays B — > D°K~ , 
D + K~, D*°K~ and D* + K~ that are well separated from the more abundant Cabibbo 
favored B — > D^n~ processes. We report values for the ratios of branching fractions R = 
B{B — > D^K~)/B(B — > D^ix~) that agree, within errors, with naive model calculation. 

We wish to thank the KEKB accelerator group for the excellent operation of the accelera- 
tor. We acknowledge support from the Ministry of Education, Science, Sports and Culture of 
Japan and the Japan Society for the Promotion of Science; the Australian Research Council 
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contract No.2P03B 17017; the Ministry of Science and Technology of Russian Federation; 
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TABLES 

TABLE I. The fit results for the numbers of D^tt~ and D^K~ events, feed -across from 
£)(*) 7r to D^K enriched sample and the statistical significance of D^K signal. 









d^k- 




stat. 






events 


events 


feed- across 


sig. 


B~ 


^D°h- 


2402.8 ± 97.8 


138.4 ± 15.5 


52.0 ± 11.4 


11.7 


B° 


-»• D+h- 


681.9 ±32.1 


33.7 ±7.3 


10.4 ±4.9 


6.1 


B~ 


-»• £>*°/r 


584.8 ± 32.4 


32.8 ±7.8 


6.8 ±4.9 


5.8 


B° 




640.9 ± 30.8 


36.0 ±7.1 


21.0 ±5.7 


7.6 



TABLE II. The resulting branching fraction ratio measurements. The first error is statistical 
and the second is systematic. 



B{B~ — > D°K~)/B(B~ — > D°ir~) = 0.079 ± 0.009 ± 0.006 

B(B° -» D+K-)/B(B° -> D+ir-) = 0.068 ± 0.015 ± 0.007 

B(B~ -> D*°K-)/B(B- -»• D*°7r-) = 0.078 ±0.019 ±0.009 

B(B° -> D*+K )/B(B° -> D* + tt~) = 0.074 ± 0.015 ± 0.006 
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FIGURES 
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FIG. 1. The AE distributions for the (a) B" -» D tt~ and (b)B° -» D+tt" samples, and 
the (c) 5~ — > D°K~ and (d)i?° — > D + K~ enriched samples, where in each case a pion mass is 
assigned to the ir~ /K~ track candidate. The points with error bars are the data, the curves show 
the results of fits that are described in the text. The open histograms are the sums of background 
functions scaled to fit the data and the hatched histogram indicates the continuum component of 
the background. 
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FIG. 2. The AE distributions for the (a) B~ L>*°vr" and (h)B° -> Z?* + vr" samples, and 
the (c) — > D*°K~ and (d)£?° — » D* + K~ enriched samples, where in each case a pion mass 
is assigned to the tt~~ fK~ track candidate. The curves and histograms are the same as those in 
Fig. |. 
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